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Abstract 

We investigate the Kosterlitz-Thouless transition for hexatic order on a free 

fluctuating membrane and derive both a Coulomb gas and a sine-Gordon 

Hamiltonian to describe it. The Coulomb-gas Hamiltonian includes charge 

densities arising from disclinations and from Gaussian curvature. There is an 

interaction coupling the difference between these two densities, whose strength 

is determined by the hexatic rigidity, and an interaction coupling Gaussian 

curvature densities arising from the Liouville Hamiltonian resulting from the 

imposition of a covariant cutoff. In the sine-Gordon Hamiltonian, there is 

a linear coupling between a scalar field and the Gaussian curvature. We 

discuss gauge-invariant correlation function for hexatic order and the dielectric 

constant of the Coulomb gas. We also derive renormalization group recursion 

relations that predict a transition with decreasing bending rigidity k. 
PACS numbers: 05.70.Jk, 68.10.-m, 87.22.Bt 
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I. INTRODUCTION 



Bilayer fluid membranes [|I],0 spontaneously self assemble when aliphatic molecules are 
dissolved in water at a sufficiently high concentration. At high temperature, these mem- 
branes have no internal order and can be modeled as fluctuating structureless surfaces char- 
acterized by a bare bending rigidity n. The bending rigidity is length-scale dependent and 
becomes zero at the persistence length = ae'^'^'^^^'^ where a is a molecular length and T 
is the temperature. At length scales less than of order C,p, the membrane is flat; at longer 
length scales, it is crumpled. 

A flat rigid membrane can have quasi-long-range (QLR) hexatic order [^j at low tem- 
perature and undergo a Kosterlitz-Thouless (KT) disclination unbinding transition to 
a disordered high-temperature phase. A fluctuating membrane can also have QLR hexatic 
order 0. Hexatic order stiffens the bending rigidity so that, rather than scaling to zero at 
long length scales, it approaches a constant times the hexatic rigidity K @J^. The hexatic 
membrane is thus more rigid than a fluid membrane, and it is said to be "crinkled" rather 
than crumpled. A fluctuating hexatic membrane can undergo a KT transition from the 
crinkled to the crumpled state. Reference discussed two possible mechanisms for the 
crinkled-to-crumpled transition: disclination melting and crumpling. The latter mechanism 
is analogous to that producing the flat-to-crumpled transition in tethered membranes |TD 
and is argued to be associated with the buckling instability of a membrane with a single 
disclination. Figure |l| shows a schematic phase-flow diagram in the {l3K)~^-{f3K,)^^ plane 
(/? = 1/T) for a 2d membrane embedded in three dimensions adapted from Ref. |^. The 
vertical line at (jSK)^^ = -njll is the Kosterlitz-Thouless disclination unbinding line of a 
flat membrane. The curved line joining the vertical line at (/Jk)^^ = vr/ll is an estimate 
of the crumpling transition obtained by equating the energy of a single positive disclination 
in a buckled membrane to its entropy. Thus, the crumpling transition in this estimate is a 
Kosterlitz-Thouless transition in a buckled membrane. This schematic phase diagram de- 
scribes qualitatively features that are in agreement with simple physical reasoning: for large 
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K, there should be a dischnation melting to the crumpled phase as temperature is increased, 
and at fixed there should be a transition to the crumpled phase as /t is decreased. It, 
however, has features that are either unexpected or unexplained. The discontinuous change 
in curvature where the melting and crumpling lines join is surely an artifact. It suggests 
that the physics of the crinkled-to-crumpled transition produced by decreasing K at fixed 
K and by decreasing k at fixed K are totally different, the latter being associated with the 
buckling instability of a membrane with a single disclination. It leaves unanswered whether 
the buckling instability line for the zero temperature membrane has any significance for a 
membrane in thermal equilibrium, which is allowed to choose disclination configurations to 
minimize its free energy and thereby to reject highly energetic configurations with an excess 
of positive or negative disclinations. 

In this paper, we present a detailed analysis of the low-temperature crinkled-to-crumpled 
transition in hexatic membranes. Our approach treats both disclinations and Gaussian 
curvature in the same real-space renormalization procedure and allows us to obtain recursion 
relations for k, K and the disclination fugacity y. Previous treatments used momentum space 
renormalization procedures to calculate the recursion relation for K and did not actually 
provide a complete set of recursion relations for k, K and y nor a prescription for doing 
so. Our procedure shows that thermally induced shape fiuctuations cause a ^-dependent 
reduction in missed in previous calculation, that leads to the phase-fiow diagram shown 
in Fig. 0. The vertical line in Fig. |l| is now curved, and the mechanism for the crinkled-to- 
crumpled phase transition in the vicinity of P, the termination of the crinkled line, is the 
same for decreasing both k and K. Our renormalization equations allow us to study the 
fiuid phase in the vicinity of the transition point P and to show that near P the persistance 
length is ^kt^'^'^'^^^^ rather than ae'^'"^/^'^ , where ^kt = aexp{b/\T - TatI^^^) is the KT 
correlation length (with b a constant), and k is the bending rigidity at length scale C,kt- 

Our calculations are strictly speaking restricted to < 1, {PK)~^ < 1 and 

{(3K)/{(3K,y = TK/k^ < 1, i.e., to the region below the curve OA in Fig. ^. We can- 
not, therefore, make any definitive statement about the interesting K oo limit, which 
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should be related to the crumpling of tethered membranes. However, we believe that the 
phase diagram shown in Fig. ^ makes physical sense beyond the region of validity of our 
calculations. In this scenario, the transition from the crinkled to the crumpled phase would 
be controlled by the fixed point P in Fig. ^for all > and n"^ > 0. The transition at 
K = oo would, however, be controlled by another fixed point. 

We begin in Sec. II with a review of how to describe tangent plane order on fluctuating 
surfaces. We pay particular attention to correlation functions of orientational order. In 
order to compare tangent plane vectors at two different points on the surface, it is necessary 
to parallel transport one of the vectors along some path F to the position of the other 
vector. When Gaussian curvature is nonzero, the direction of a parallel transported vector 
depends on F even when there are no disclinations present. Physical correlation functions 
are invariant with respect to local coordinate transformations and for a particular membrane 
shape and distribution of disclinations depend on F. When correlation functions are averaged 
over shape and position of disclinations, the dependence on F vanishes. 

In Sec. Ill, we discuss various models for hexatic membranes. We begin with the Hamil- 
tonian for hexatic membranes expressed in terms of the orientational angle 6 and a height 
variable. We then transform this model into a Coulomb gas model in which "charge" density 
arises both from disclinations and from Gaussian curvature. There are two Coulomb-like 
terms in this Hamiltonian: one which is zero when the local disclinations density equals the 
local Gaussian curvature and one, arising from the imposition of a covariant cutoff |12||, that 
couples Gaussian curvature to Gaussian curvature. Finally, we transform the Coulomb-gas 
Hamiltonian to a sine-Gordon Hamiltonian with a term coupling the sine-Gordon field to 
the Gaussian curvature with an imaginary coefficient. The latter term is analogous to the 
dilaton coupling of string theory ||14|| . 

In Sec. IV, we relate the dielectric constant and hexatic rigidity to correlation functions 
of the disclination- Gaussian curvature density. We show in particular that the renormalized 
hexatic rigidity appearing in orientational correlation function is the same as that calculated 
from the free energy. We also calculate the charge-density correlation functions. 
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Finally in Sec.V, we derive renormalization group recursion relations for K, k, and 
the disclination fugacity y. These equations show that height fluctuations renormalize the 
hexatic rigidity in the absence of disclinations. This renormalized rigidity then flows under 
renormalization in the presence of disclination in exactly the same way as the rigidity of a 
flat membrane. If the initial height renormalized rigidity is less than the critical rigidity for 
a flat membrane, there is no rigid phase. Thus, height fluctuations can destroy the crinkled 
phase. 

In Sec. VI, we review results and discuss some unanswered questions. In particular, we 
address the effect of asymmetry in the energies of positive and negative disclinations [l^ on 
our results. 

The focus of this paper is on the nature of tangent-plane order and the Kosterlitz- 
Thouless transition on fluctuating surfaces. It will not present a complete derivation of the 
RG equations used in Sec.V because they involve some two loop graphs, which should be 
treated in a sophisticated regularization procedure. In a companion paper, we will derive the 
complete recursion relations from the sine-Gordon model using a real space regularization 
procedure generalized from that used by Amit et al. W^ to treat the flat space problem. 



II. DIFFERENTIAL GEOMETRY AND TANGENT-PLANE ORDER 

We are concerned with the nature of tangent-plane order on fluctuating membranes. In 
this section, we will review relevant concepts in differential geometry, mostly to establish 
notation. We will also discuss how to describe long-range order (or lack thereof) on a metric 
with fluctuating curvature. 



A. Differential Geometry of a Plane 

Points [p!^^-[T9| on a two-dimensional surface embedded in three-dimensional Euclidean 
space are specified by a three-dimensional vector R(u) with components Ri{u), i = 1,2,3, 
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as a function of a two-dimensional parameter u = (w^, u^). Covariant tangent- plane vectors 
are then defined as 

ta^daR, a = 1,2, (2.1) 

where da — d/du°'. We will use Greek letters a, /3,7, ... to denote components of covariant 
and contravariant tangent-plane tensors and Roman letters k, ... to denote components 
of vectors and tensors in Euclidean space. The metric tensor is 

gafS = ta • tp. (2.2) 

Its inverse gf"'^ satisfying 

g'^^gp, = (2.3) 

allows us to define contravariant tangent- plane vectors t° = g^'^tp satisfying t" ■ t/3 = 5 p. 
Any vector V in the tangent plane can be expressed as V = V^ta = Vat"' where = ta • V 
and V°' = t" • V = g'^^Vp are, respectively, the covariant and contravariant components of 
V. A unit normal N to the surface can be constructed from ti and t2 : 

N = (2.4) 

The curvature tensor is then 

= N • dadp-R. (2.5) 

Prom the curvature tensor, one can construct the mean curvature, 

= + (2.6) 

2 2 2 Vi?i i?2/ ^ ' 

and the Gaussian curvature, 

5 = deti^^ = i--i-, (2.7) 

where Ri and R2 are the principal radii of curvature at the point of the surface in question. 
The integral of the Gaussian curvature is a topological invariant, 
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I d^u^S = 47r(l - 77) = 27rx, 



(2.8) 



where r] is the number of handles and x — 2(1 — rj) is the Euler characteristic. In the Monge 
gauge, u = {x,y) and R(u) = (u, /i(u)), and the metric tensor g^jj is written as 



gap — ■ dfjK — 
and the curvature tensor K^p is 

Kap = ^- D^D^R 



1 + {dxhY dxhdyh 
^ d^hdyh 1 + {d^hfv ^ 



^ d^dxh dxdyh^ 



(2.9) 



VTT(V^ 



y dyd^h dydyh J 



(2.10) 



where {Vhf = {d^hf + {dyhf. 

The anti-symmetric tensor ^^13 will be particularly useful in what follows. It is defined 

via 

7a/3 = N ■ (t« X tp) 

gai9i32 — ga2gfil 



|ti X 



(2.11) 



where g — det g^p and is the anti-symmetric tensor with 612 = —^21 — 1- The contravari- 

ant tensor 



= N • (t" X t'^) 

equals eapj ^ and satisfies 7"'^7/3a' = ~^a'- Finally the mixed tensor 



(2.12) 



1 = g la'n 



(2.13) 



rotates vectors by 7r/2 since Kt"/?^^ = 7a/3^"^^ = and T^^T^^Ta^V^?/ = 

We are interested primarily in order in the tangent plane of unit (or fixed) magnitude. 
For this, it is useful to introduce orthonormal tangent-plane basis vectors ei and 62 satisfying 



ea • 66 = boh-, N • ea = 0. 



(2.14) 



A tangent vector V can be expressed in the basis {ei, 62} as well as that defined by the co- 
variant or contravariant vectors : V = V^e^ where Va — e^- V. We will use Roman subscripts 
a,b,c, ... to denote vector and tensor components with respect to the local orthonormal ba- 
sis. Covariant derivatives are derivatives projected into the tangent plane. Components of 
the covariant derivative of a vector V relative to the orthonormal basis are 

DaVa = e„ • (9a V) = + ■ ^aCfeH 

= daVa + eabA^Vb, (2.15) 

where 

Aa^ei-dae2 (2-16) 
is the spin-connection whose curl is the Gaussian curvature : 

r^d^A^ = S. (2.17) 
We will also find it useful to use a circular basis defined by the vectors 



e± = -^(ei±ie2)=e;, (2.18) 



satisfying Sa ■ el — Sab with a,b — ±. In this basis, V — 14e*, and the covariant derivative, 

DaV± = e± • daV = daV± + e± • daelVa 

= {daTtAa)V±, (2.19) 

has a particularly simple form. 

B. Vector and Tensor Order 

A vector order parameter S that is restricted to lie in the tangent plane of a surface can 
be written as S'^ta or Sa^a- If S is a unit length vector, it is conveniently expressed in terms 
of its angle in the local orthonormal basis as 



S = COS 6161 + sin 9e2 = SaBa, (2.20) 

where Si — cos 9 and 5*2 = sin^. The unit vector S can thus be written in the circular basis 
as 

S = V2Re* (2.21) 

where 

* = Ve- = e^^e_ (2.22) 

is a complex order parameter with unit length ^ • ^* = 1. 

We now turn to tensor tangent-plane order. The simplest nontrivial tensor order param- 
eter is a symmetric-traceless tensor Q with Cartesian components 

Qij — Qab^ai^bj ■ (2.23) 

The traceless constraint implies Qa — Qaa — 0. Q is a uniaxial tensor with a principal axis 
lying along a unit vector N = riaGa in the tangent plane. If we require that Q have a fixed 
magnitude defined by TrQ^ = 1, then we can write 

Qab = V2{nanb - ^Sab) (2.24) 

or 

1 

Qij = —^[cos29{eueij - 62162^) + sin 26'(eiie2j + 62161^)]. (2.25) 

This tensor, like the vector S, can be expressed in terms of the real part of a complex tensor. 
Introduce the direct product tensor 

*2 = * <2) * = -026- (8) e_ = e^^'^e^ ® e_ (2.26) 

with components — e^*^e_je_j. Then 

Q = V2Re*2- (2.27) 



Generalization of this construction to higher order tensors is straightforward. Let 

*p = * . . . (g) = ^pe_ (g) • • • (8) e_ (2.28) 

with %l)p — e^'^ be a pth rank complex tangent-plane tensor. If we define the inner product 
of two tensors by 

ySi . S>* = \[i S>* \I/ . \I>* = IvEr |2 — 1 (0 2Ql 

^^l■■■^p^^-^■■■^pJ ^p \^p\ y^.^nj 

then a real tensor, symmetric under interchanges of all indices and traceless with respect to 
all pairs of indices, is 

Qp = V2Re *p. (2.30) 

The tensor = ^' describes vector order such as is present in the smectic-C phase in 
which the long axes of molecules comprising the membrane tilt relative to the normal N. 
^2 describes tangent-plane nematic order with inversion symmetry, describes hexatic 
order. ^4 would describe "4-atic" order, etc. Thus for general p describes what we call 
"p-atic" order. 

C. Tangent-plane Correlations and Parallel Transport 

In flat space, the basis vectors 61^2 or e± are independent of spatial coordinate. Informa- 
tion about the existence of long-range order and about p-atic order parameter correlations 
is contained in the correlation function (^p(x) ■ 'J'*(x')) = (e"*^^^^''')"^^^))). Basis vectors 
at different points on curved surfaces are not identical (or even in the same plane) and the 
simple dot product of vectors at different points does not carry information about tangent- 
plane order. In order to compare order parameters *p(u) and **(u') at two different points 
u and u', we need to parallel transport the order parameter at u along some path T to u'. 
The parallel transported tensor *p(u, u') is now in the tangent plane at u' like *p(u') and 
we can take its dot product with **(u'). Thus the correlation function. 
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G,(u, u') = (*,(u) : *;(u')) ^ (*,^(u, u') • *;(u')), (2.31) 

where the ":" means parallel transported inner product, is what is appropriate for describing 
correlations in tangent-plane order. Note that this correlation function depends on the path 
r when the membrane has nonzero Gaussian curvature. We will find, however, that it 
becomes independent of F on nearly flat membranes after averaging over height fluctuations. 

The parallel transported order parameter can be expressed in terms of the spin- 
connection. Let V(u) = Va(u)ea(u) be a tangent vector at u and let V|l(u + 6u) = 
Vj{u + Su)ea{u. + Su) be the vector V(u) parallel transported to a nearby point u -|- Su. 
Then by definition yjl(u + 5u) = K(u) + SVa{u) with 

6Va{u) = -e„5A„H(u)5ii". (2.32) 

Alternatively in terms of the circular basis, 

6V±{u) = ±iV±{u)Aa{u)6u'', (2.33) 

or 

VX(u, u') = e±'/u'^«(")''«"-(4(u). (2.34) 
Applying this result to the vector * = V'^-; we obtain 

*^(u, u') = e'/u' ^«(")<i""*(u), (2.35) 

or 

*^(u, u') = e*/u'^"(")'^-"e''(")€_(u). (2.36) 

Thus we have 

Gi(u, u') = (^e-'^o(u'yeiuyj:' A^iu)du'')^ ^2.37) 

This function is invariant under changes in the local coordinate system (i.e. under nonlocal 
rotations of the vectors ei and 62). The generahzation to p-atic order is straight forward 
since by construction *p(u) : **(u') = (*(u) : **(u'))^' so that 
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Again, this is invariant under coordinate system changes. 



(2.38) 



III. MODEL HAMILTONIANS 

In this section, we will derive various equivalent representations for the Hamiltonian and 
associated partition function describing p-atic order on a fluctuating surface. 



A. Fluid Membranes 

It is now well established that the long wavelength properties of a fluid membrane are 
well-described by the Helfrich-Canham Hamiltonian TYhc |2C,21,1|], which can be expressed 
as a sum of three terms, 

Hhc = 'Hk + Hg + Ti.a- (3.1) 
The first term is the mean curvature energy, 

Vh V 



(3.2) 



2 J l + (V/i)V 

where H = K^, and the second form is valid for the Monge gauge. The second term is the 
Gaussian curvature energy, 

= ^i^G J d^Uy/gS, (3.3) 

where S = deti^^. This term is a topological invariant depending only on the genus of the 
surfaces. Since we will consider surfaces of fixed genus, we will drop this term. Finally 

H, = a j d\^ (3.4) 

is the surface tension energy. We are mostly interested in free membranes for which the 
renormalized surface tension obtained by differentiating the total free energy JF with respect 
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to the total surface area A (o"_r = dJ^/dA), is zero. Since there are entropic contributions 
to an as well as contributions vfrom internal order, the value of the bare surface tension 
a will have to be adjusted to keep aR zero. In what follows, we will ignore Ha with the 
understanding that it is really present if we want to keep track of how actually becomes 
zero. 

The partition function for a fluid membrane, 

Zfi = J r'R(u)e-^^H° (3.5) 

is obtained by integrating over all physically distinct realizations of the surface, which is 
specified by the vector R(u). This means we have to specify a gauge or parametrization 
for R(u). Thus the measure PR(u) for distinct surface configurations contains a Fadeev- 



Popov determinant p2[. In addition, it in general contains a factor to correct for the fact 
that different surface configurations arising from a fixed parametrization surface can have 
different areas 



B. p-atic Order 

As discussed in the preceding section, the complex tangent-plane tensor distinguishes 
a fluid membrane from one with p-atic tangent-plane order. Since has a fixed magnitude 
and there are no external fields aligning along a particular direction, the lowest nontrivial 
contribution to the energy associated with arises from its gradients. 

He = \Kp J d'u^g"''D^% ■ (3.6) 



where Da is a covariant derivative. Using Eq. ( p.l9|) and Eq. (|2.28| ), it is straight forward 
to show 

Da^p = ip{dJ-Aa)^p (3.7) 

and 



13 



'He = \K j d\^g'^l'{dj - A^){dp9 - Ap), (3.8) 



where we set K = Kpp^. Eq. ( |3.8| ) is the simplest energy arising from p-atic order. It is 
rotationally isotropic, and it describes correctly the lowest order gradient contribution to 
the energy for all p > 3. For p = 1 or p = 2, however, this energy should have anisotropic 
contributions [^,^. In this paper, we will ignore these anisotropics. 



The p-atic order parameter can have disclinations of strength q = 2'K(k/p) where k is 
an integer |2^. A disclination at u = u,^ with strength gives rise to a singular contribution 



9^^°^ to satisfying 

j( du^dj''''^ = q„, (3.9) 
where F is a contour enclosing Uj^. Thus, in general daO = dad' + Va where 6' is nonsingular, 



= dj''"^^ and 



7°^9,t;^ = n(u), (3.10) 



where 



n(u) = ^ Vg,.5(u - u^) (3.11) 

is the disclination density. The vector f „ can always be chosen so that it is purely transverse, 
so that -Daf " = 0. In the p-atic Hamiltonian, d^fi always occurs in the combination daQ — A^- 
The spin-connection A^ can and will in general have both a longitudinal and a transverse 
component. However, one can always redefine Q' to include the longitudinal part of A^- This 
amounts to choosing locally rotated orthonormal vectors ei(u) and 62 (u) so that -D„A° = 0. 
Thus we may take both f „ and A^ to be transverse and the p-atic Hamiltonian, 

ne = ^Kj d\^g^^{dj + Va- Aa){dpe' + vp- Ap) (3.12) 

= + 7i_L, 

can be decomposed into a regular longitudinal part, 
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Hji = ^kJ d\^g^^dj'dpe', (3.13) 

and a transverse part, 

n^ = ^Kj (fu^g'^^iv^ - A^){vp - Ap), (3.14) 

where it is understood that -Da(f " — A") = 0. 

It costs an energy edk) to create the core of a dischnation of strength k. (We assume 
for the moment that the core energies of the positive and negative dischnations are the 
same. See, however, Ref. |]T3| and the summary section.) Thus, partition sums should be 
weighted by a factor yu = e"^^''^^-' for each dischnation of strength k. Since ec{k) ~ k'^, we 
may at low temperature restrict our attention to configurations in which only configurations 
of strength ±1 appear. Let N± be the number of dischnations of strength ±1 and let Ui,± 
be the coordinate of the core of the dischnation with strength ±1 labeled by u. The p-atic 
membrane partition function can then be written as 

Z{K,K,y) = Ti^y^ j DR j m-^^'^e"^^*, (3.15) 

where y = yi, and = + A^_. Tig depends on all of the dischnation coordinates U;,± 
where = 1,2, ■ ■ ■ , N±, and 



Tr. = E V+7v_,pxTni^n / / (3-16) 

where a is a molecular length. The Kronecker factor 5Ar_^+7v_,px in Tr^ imposes the topological 



constraint that the total dischnation strength on a surface with Euler characteristic x 
be equal to x- Thus, on a sphere with x = 2, + = 2p. If p = 6 for example, there 
must be 12 more positive than negative dischnations. On a nearly fiat surface, the Coulomb 
interaction effectively restricts to equal N_, and we use Eq. (|3.16|) with x = even 
though X for an open surface is strictly speaking equal to one . 



C. The Coulomb Gas Model 



The p-atic model of Eq. (|3.15|) can easily be converted to a Coulomb gas model using 



15 



r''d^{vp-Ap)=n-S = p, (3.17) 



which follows from Eq. ( p.l7|) and Eq. ( |3.10|) . The quantity p = n — S is a, "charge" density 



with contributions arising both from disclinations and Gaussian curvature. Equation (|3.17|) 
implies 

v^- A^ = -^jDp^p, (3.18) 

where = D'^Da = {^/ y/g)da^/gg°'^df^ is the Laplacian on a surface with metric tensor 
ga/3 acting on a scalar. Recall [Eq. ( p.l3|) ] that 7q,^ rotates a vector by 7r/2 so that Va — 



is perpendicular to Df^{—Ag) ^p and is thus manifestly transverse. Using Eq. ( p. 18 ) in 
Eq. ( |3.14|) , we obtain 



Z = Tr.y^l PR I p^'e-^^--'3^^r^^^ (3.19) 

where 

'Hc = \k j dVu'^pin) (-^] v^p(u') (3.20) 

\ 9 / uu' 

is the Coulomb Hamiltonian associated with the charge p. The longitudinal variable 6' 
appears only quadratically in Tiy and the trace over 6' can be done directly: 

J D^'e^^^ii = e-'3^^ (3.21) 



where 



= 8^ / '^'''^ " i / d'^d\'^S{u) (-^) ^S{u') (3.22) 

\ 9 / uu' 



is the Liouville action [lT2| , p!9| arising from the conformal anomaly. The first term in this 
expression is proportional to the surface area and can be incorporated into the surface tension 
energy by shifting the surface tension a. We assume, as discussed earlier, that the total 
surface tension is zero. We will, therefore, ignore this term in what follows. The second 
term can be viewed as a Coulomb interaction for Gaussian curvature but with a negative 
coupling constant, i.e., an imaginary charge. 
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The Coulomb gas partition function can thus be written 



where 



Z = Ti^y^ J VRe~ 

= Tr.y^ J PRe-/3W-^^cT ^ (3 23) 

\ 9 / uu' 

+ J d\d\'^n{u) ^n{n') 

\ 9 / uu' 

-K J d\d\'^n{u) (-^) ^/g'SiW) 

\ 9 / uu' 

= 'Hss + 'Hnn + Ti-nS (3.24) 

is the "total" Coulomb Hamiltonian in which there are two distinct charge densities, p and 
S, or, equivalently, n and S. The coupling constant for Gaussian curvature interactions in 
the absence of disclinations is 

A" = A--^. (3.25) 

This is the shifted coupling arising from the Liouville term discussed in Ref . and . The 
coupling constant for disclinations and between disclinations and Gaussian curvature are 
not shifted by the Liouville term and remain equal to K. 

D. The sine-Gordon Model 

The Coulomb gas model can be converted following standard procedures into a sine- 
Gordon model. The first step is to carry out a Hubbard-Stratonovich transformation on 
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g-/3Hc = g/9HL J p^g-l/CW) / d^u^d'-<t,dc,'t>^i J cfu^p<t>^ (3.26) 

where the Liouville factor e'^^^ is needed to ensure that e~^^'^ be one when p = 0. Inserting 
this in Eq. ( |3.23| ), we obtain 

Z = Ti^y^ J VRV^e'^'^'^e-'^'^^e' / '^'"^/9(«-^)<^, (3.27) 

where 

1 r,2. 



/9H<^ = ^ y d'u^d''(l)d^(j). (3.28) 

The only dependence on dischnations is now in the term hnear in n. Thus to carry out Try, 
we need only to evaluate 



(/^^^.-.<w/.)" 

V ^ f ^e»^PX /^^ /■ ^!V^gi{2^[<^(u)/p]-..}\ ^ (y f ^!V^g-i{27r[0(u)/p]-<..} 

A^^!A^_! i 27r \ i j \ J 



N+,N- 

dui 



g«^PXg(2j//a2) J a!2„^cos[27r(0/p)-i^] 

27r ■ ^ ' ' 



Thus 

duj 



Z = j — J V(f) y"DRe-^''^«e"'^^'*e*'^P^e(^^/"''') / d2«^cos[27r(0/p)-Hg-i / <fu^Sci,_ (3.30) 

We can now change variables, letting (p = (p/27r)(0' + u). The term linear in the Gaussian 
curvature then becomes 

- i J d^Uy/gS^{u + (j)') = -iujpx " / d'^u^/gScp' , (3.31) 



where we used Eq. ( ^.8| ) relating x to the integral of the Gaussian curvature. Thus the 
cj-dependent part of this term exactly cancels the iujpx term arising from the integral rep- 
resentation of the Kronecker delta. The integral over a; in Eq. ( [j.3(J| ) is now trivial, and we 
obtain 



Z = y y PRe^^^^e-^, (3.32) 

where 

^=^{£}' J d^^g^'^d^cl^dpcl) - J y Su^cosc^- J d\^Scl> (3.33) 

is the sine-Gordon action on a fluctuating surface of arbitrary genus. The first two terms 
of this action are the gradient and cosine energies present in a flat space. The final term 
provides the principal coupling between and fluctuations in the metric. It is analogous 
to the dilaton coupling of string theory though here the coupling constant is imaginary 
rather than real. Note that the Liouville action is not explicitly present in Eq. ( p.32|) . The 



requirement that the cutoff be applied in a covariant fashion is, however, still present. 
Thus, if (7 = 0, the integral over will yield the Liouville factor. When y is not zero, care 
must be taken to implement any cutoffs in integrals in a covariant fashion. 

IV. THE DIELECTRIC CONSTANT AND HEXATIC RIGIDITY 

In the preceding section, we derived several equivalent expressions for the partition func- 
tion of a fluctuating membrane with internal p-atic order. In this section, we will derive 
an expression for the inverse dielectric constant associated with the charge density p and 
show that it is equivalent to the renormalized p-atic rigidity controlling the gauge-invariant 
p-atic correlation function of Eq. ( |2.31| ). We will begin by reviewing the derivation of the 
longitudinal dielectric constant in flat space. We will then derive the dielectric constant and 
renormalized rigidities on fluctuating surfaces. 



A. Flat Space Dielectric Constant 

Consider a system with an internal charge density p and fixed or controllable external 
charge density pe- To be concrete, pe could be the free charge at the metal electrodes of a 
capacitor. The total charge density is = P + Pe- The electric field £ is determined by the 
total charge: 
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V • = apT- (4.1) 

Here a sets the units of charge (for example a — in mks units). The displacement vector 
T> — ea'^E, where e is the dielectric constant, is controlled by the external charge: 

V • P = V • P|| = pe, (4.2) 

where T>\\^ is the longitudinal part of T>. The electric potential can be introduced in the usual 
way: £ — — V0t, = — a;~-^V0e and 0t = + 0c where — = ap and — V^0e = ctpe- 

The inverse longitudinal dielectric constant is obtained by expanding the free energy to 
second order in T>\\^ — — V0e- 

AJT = 1q;-i j d'^xd'^x'e^\x, x') V0e(x) • V0e(x'). (4.3) 

A wave-number dependent dielectric constant can be obtained if e|j"^(x, x') is translationally 
invariant or by carrying the usual Maxwell averaging procedure. The result is 

where V — J dfix is the volume of the system. We can use this formula to relate e|[''^(q) to 
correlations of the charge density. The Coulomb Hamiltonian can be written as 

= j d''x{V(f),f + d'^xpcf) + j d'^xpcf),. (4.5) 

Expanding — —Tlne~^^^°~^''^'^\ where Tio is the non-Coulombic contribution to the 
Hamiltonian to second order in V0e) we obtain 

A.F = ^a-' J d'x{W<Pe)' -^pj d'^x0e(x)(/>e(x')(p(x)p(x')) (4.6) 



and 



where 
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Cpp(q) = IpI d''xd''x'e-'^<^-^\pi^)pix')) (4.8) 

is the density correlation function. Eq. ( [4.7[ ) is the well-known expression relating e^^ to 
the charge susceptibility Xpp = PCpp. In what follows, we will use the notation introduced 
in Eq. (|4.8| ) for arbitrary variables, 

Cab{cD = / dWx'e-^''-(^-^')(A(x)i?(x')), (4.9) 

where we have set = 2 and V = Ab- 



B. Dielectric Constant on a Fluctuating Membrane 

We can now determine the dielectric constant associated with the charge p = — 5* on a 
fluctuating membrane following exactly the procedures outlined in the preceding subsection. 
We impose external charges pc to create a constant slowly varying external potential 0e- 
These charges lead to a total charge density px = p + Pe- They should not, however, 
change the Gaussian curvature because otherwise they would change the Liouville energy. 
Thus Pe must arise from disclinations which we can, for concreteness, place at the boundary 
of the membrane like the charge on capacitor electrodes. In analogy with flat space, we 
can introduce external and induced electric potentials 0e and satisfy — Ag0 = p and 
— Ag0e = Pe with 0T = + 0e the total potential. The Coulomb Hamiltonian can then be 
written as 

+ \j d\y/^p(l) + j d\^p(t)e. (4.10) 

We are interested in the long wavelength dielectric constant for our presumed macro- 
scopically flat (or nearly flat) membrane. For this purpose it is most convenient to use the 
Monge gauge in which u = x = (x, ?/) measures position in the average plane of the mem- 
brane. In this case 0e(u) 0c(x) is a slowly varying function of x, which, in the limit of 
fixed charge density at membrane boundaries, grows linearly with x. The long wavelength 
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dielectric constant for our nearly flat membrane can be obtained using Eq. ( [4 .41) , where q is 
the Fourier variable dual to position x, and Eq. ( [4 .91 ) for TIq. Using Eq. (f4.9| ) and expanding 
JF to second order in 0e, we obtain 

I d\ I Ci\'(v/^p(x)y^p(x'))0e(x)0e(x') (4.11) 

and 



^r^(q) = ^(v^^"')-^c^pp(q) (4-12) 

is the charge density correlation function for 

P = V9P (4.13) 

rather than p. p is in fact the effective charge density for the coarse grained flat surface 
(with metric tensor equal to unity) because the total charge in an area patch d'^x should 
not change under coarse graining. The quantity {^g'^^) is equal to 5°^ as can be verified 
easily to lowest order in T in the Monge gauge or more generally using the conformal gauge 
Thus we have 

h\^) = ^-^CM- (4.14) 

Cpp(q) is the charge density correlation function evaluated in the presence of all interactions, 
including the Coulomb like TYl. We will return at the end of this section to its evaluation. 

C. Hexatic Rigidity 

In this section we will show that the renormalized hexatic rigidity Kji is equal to ae|j"^. 
To calculate Kr, we introduce a shift uJa in daO' via da_0' — > daO' + in Eq. ( |3.13| ) and 
evaluate the free energy to second order in oja- Using Eq. ( [4.11| ), we obtain 
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^j3K^ I d\ f rf\'(v^K - - A'^)u;^u;p). (4.15) 



2 



As discussed previously, 6' decouples from other variables. For any given realization 
of the surface {06') = {T / K){—A~^). We can therefore combine the first two terms in 
Eq. ( [4.151) to obtain 



-\PK^ j d\ J d\'{^v^uj^^v^ujp), (4.16) 



where 

. . _ ^ (4.17) 

is the transverse projection operator and 

(4.18) 



is manifestly transverse from Eq. ( p.l^ ). Thus, as expected, AjF depends only on the 
transverse part of oja- 

A longitudinal function such as daipe maintains its form under coarse graining. A trans- 
verse function does not because the direction of a transverse function depends on the local 
value of the rotation matrix 7°^. Thus to coarse-grain our transverse function C(Jq,(x), we set 
it equal to 7a^9/30e- In the Monge gauge c?/30e has the same value regardless of the particular 
configuration of the membrane. Then using 7"q,/7°'^ = —^'p, we obtain 

-\PK^ j d\ J rf\'(v/^p(x)v9p(x))</>e(x)(/)e(xO. (4.19) 
Setting ea/3di3(j)e = ^^3(^0/3 — dadp/d"^) where is the coarse grained cUq, we find 
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^ Ab du^dujl 

= Ke^^\ (4.20) 

where Ab is the base area and e\\ = e||(g = 0). Because / (fxp = 0, Cpp(q) tends to zero as 
g — i> in the ordered phase. We can, therefore, define 

= liml(C^^(q)-C^^(0)) 

= li^o^i^ / d'xd'x'ie^^-^-^-'^ - l)(p(x)p(xO) 

= -i/rf2x|xHp(x)p(0)), (4.21) 

where in the final step we used the fact that (p(x)p(x')) depends only on x — x' in transla- 
tionally invariant systems. Thus 

Kr = K- PK^B. (4.22) 

Cpp(q) is derived in Eq. ( ^^) . 

D. The p-atic correlation function 

We have argued that the renormalized rigidity is determined by the dielectric constant 
associated with p rather than that associated with some other combination of the inde- 
pendent charge densities n and S. We will now show that it is indeed this rigidity that 
controls the long distance properties of the p-atic correlation function Gp{u, u'), [Eqs. (p.31|) 



and ( p.38| )]. Our calculations follow closely those in Ref. for correlation functions in fiat 
space. To keep things simple, we calculate only to lowest order in height fiuctuations in the 
Monge gauge. 

As discussed in Sec. II, the angle 9 in Eq. ( |2.38| ) can be broken up into a regular part and 
a singular part: ^(x) — ^(x') = ^'(x) — ^'(x') + ds"Va- In addition, we can always choose 
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a gauge so that Aa is purely transverse. In this case 6' decouples from other variables, and 
we have 

G,(x - x') = (e-i?'^«^'W-^'('^'))^)«e-^^r''^^"(^"-^'^))R,., (4.23) 

where ((^'(x) — 9'{'x')Y)g is evaluated for a fixed metric tensor gap determined by R, and 

{A)n,. = Z-^Tr,y^ J VRe'^^'^-^'^''^ A. (4.24) 

To the order we are working, we can factorize Gp(x) into a "spin- wave" part arising from 6' 
and a transverse part: 

Gp(x) = aw(x)a(x), (4.25) 

where 

G3.(x) = e-^f^((('^'W-^'(-'))^>.> (4.26) 

and 

G,(x) = e-y^^^^^^'^ = e-^p'9v(x)^ (4^27) 



where AO^ = ds'^iva — A^) and the double brackets in Eq. ( [4.26|) refers to an average 



over 9 followed by an average over R. In the Monge gauge, to the order we are working, 

((^'(x)-^'(x'))\ = ^^(x-x'), (4.28) 

where 



Thus 



Gs„(x) = |x|-P'/2./3K_ (4^3Q) 



To evaluate G'v(x), we use Eq. (|3.18|) to lowest order in the Monge gauge to obtain 
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A^v = ea/3 / d^ypiy)fiy), 



where 



Then using the Cauchy-Riemman relation, 

where ^±(x) = (l/27r) tan(?//x), we find 

/(l/) = -(^±(x'-y)-^x(x-y)) 



and 



a(x-x') = I AA'(p(y)p(y ))/(y)/(y') 



rf'2/rfV(p(y)p(y'))(/(y)-/(y'))' 

Setting y = R + r/2 and y' = R — r/2, and expanding in r, we obtain 

G'v(x-x') = -i 1 dV(p(r)p(0))r„r;3 / d'RdJ{R)dpf{K). 
Then using Eqs. {^^^, ( TO) , (H) and V|^(s - R) = 5(s - R), we find 



Gv(x - x') = 2B[g{^ - x') - 6;(0)] 



25, |x-x'| 

— m 

2ti a 



(4.31) 



(4.32) 



(4.33) 



(4.34) 



(4.35) 



(4.36) 



(4.37) 



where we have treated ^(0) as a non-divergent constant because (x — x') is restricted to be 
greater than a. 

Combining Eqs. (Q), ( ^ ), (I^SOD and (WM) , we obtain 



X 



(4.38) 
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with 



iPKn)-' = iPK)-' + B. (4.39) 
This equation agrees with Eq. ( |4.22| ) to lowest order in B. 

E. Evaluation of the Dielectric Constant 

As we discussed at the beginning of this section, the correlation function Cpp(q) must be 
evaluated in the presence of all interactions including that arising from the Liouville term. 
In this subsection, we will derive a general expression for Cpp(q), which we will evaluate in 
the low-temperature, low-fugacity limit. We will also discuss effective Coulomb interactions. 
Both K and K have units of energy. Thus, both {j3K)~^ and {j3K)~^ are unitless measures of 
temperature. In addition either of these quantities multiplied by any function of the ratio 
K/n is also a measure of temperature. In our low-temperature expansion, we will require 
Al3K/{l3Kf = A{T/k){K/k), where A smaller than l/(27r), as well as (pn)-^ and {(3K)-^ 
to be small. Thus, we will not consider the limit K oo. We will, nonetheless, be able to 
explore a large region of the phase diagram. 

The evaluation of Cpp(q) is obtained most directly in terms of diagrams. Figure |^ 
shows the basic ingredients of a diagrammatic perturbation theory: the height propagator 
Ghh{q) = {pnq'^)~\ the three Coulomb interactions [Eq. (|3l2^) ] K'/q'^, K/q^, and K/q^ 
coupling, respectively, S to S, h to n and S to n, the h — h propagator when Tins is 
zero, and non-linear vertices arising from the curvature. The total Coulomb Hamiltonian of 
Eq. |3.24| contains S — S and n — n terms and a part 



Tins = -K J d'^u^/gn (^~^^ v^'^' 



(4.40) 



coupling n to S. We will first calculate correlation functions when Tins is zero. Then we 
wiU calculate the effects of turning on Hns- If 'Hns is zero, then C^g where S = ^/gS is zero, 
and the Gaussian curvature correlation function Cgg is simply that of a Coulomb gas with 
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coupling K' . It can be represented as a sum of polarization bubbles shown in Fig. ^ linked 
by a Coulomb interaction as shown in Fig. ^. The result is 

To lowest order in /i, the Gaussian curvature is 

S = hsj^hV^h - didjhdidjh), (4.42) 

and, to lowest order in (/^k)^^, P(q) is determined by a diagram (a) in Fig. |: 

11 f d^k {qx k)4 



2{/3k)^J (27r)2A;4(k + q)4 
3 g2 



(4.43) 



327r (/3k)-2' 

Note that V{q) is proportional to even though the numerator in the integrand is propor- 
tional to g^. There is an infrared singularity in the integral over k that converts the to 
the Higher order contributions to V are shown in Figs. ^ (b)-(f). The contribution to V 
from Fig. ^ is of order q'^j3K/{j3K)^, i.e., of order (3K/{(3k)'^ smaller than Eq. (OS). The 



diagrams of Fig. ^ and ^ are smaller by another factor of (3K/ {(3k)'^. The diagram in Fig. 

which arises from the factor and nonlinear terms in 5, is of order {(3n)~^ times Eq. 
( [4.43| ). It, however, has an infrared divergence, which would have to be handled with care 



if we wished to calculate to the next order in T. The diagram in Fig. |f is of order (/?k)~^. 

When TinS = 0, the disclination correlation function Cnn{q) can be calculated as a power 
series in y. To lowest order in T, we can replace y/g by 1 in n. Then to lowest order in y, 

ClM = j^^^y'j J ^^(n(q)n(-q))e-^«'=-/'^--/^^"", (4.44) 

where n(q) = (2TT/p){e^'^'^+ — e~^'^'^^). Therefore 



p2 



y V / drr^-^^"'' + 0{y^). (4.45) 

J a 
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When TinS is turned on, there will be additional contributions to the polarization bubble 
such as shown in Fig. ||. The correlation functions Cnn, Cgg and C-^ can be obtained from 
the diagrams in Fig. |^: 

1 - {P'Kyq^)CLCl^ 
V 

(4.47) 



1 + {PK'/q^)V - p^K^VCi 



C- = ^ ^^"^^^ f4 48) 

where is given by Eq. (|4.43|) with 'P(q) corrected by contribution such as that in Fig. 



^. Combining Eqs. (|4.471 ) and ([4.48|) , we obtain the charge- density correlation function 
determining K^: 

CM = CU^) - 2Q^(q) + CsM) 

_ C9,(q) - 2{f3K/q')CU^)CssicD + Q^(q) 

1 - i/3^Kyq^)CUci)CsM 
_ ^(q) + - mK - K')/q^)V{ci)] 

1 + m/q^)V{^)[{K'/K) - (/3«:/g2)C?,(q)] " 

In the limit of low-temperature and low-fugacity, this reduces to 



(4.49) 



Cpp(q) ^ V{c,) + ClM [l - 2^P(q)) + ■ ■ • , (4.50) 

where P(q) is given by Eq. (|4.43|) and C°^(q) by Eq. (|4.45|) . The second term in the 
coefficient of C^f^ is of order PK/^PkY and will be neglected. 

We can now calculate the dielectric constant or renormalized rigidity to lowest order in 
temperature and fugacity using Eqs. (|]21|), (|0^ ), (|0^ ), (|^) and ( ^30| ). 



47r — 

(3Kr = I3K- {(3Kf—y^ / drr3-2-/^^, ^4 5^^, 

P J a 

where 



pK = pK-J-(—]\ (4.52) 



327r \ K J 
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we assume {3/327i)PK/{Pk)'^ ^ 1 so that 



(4.53) 



Thus, the equation for (3Kji is identical to that for a rigid flat membrane but with K 
replacing K. 



V. RENORMALIZATION EQUATIONS 

To determine the properties of the Kosterlitz-Thouless transition on a fluctuating mem- 
brane, we need to determine the renormalization flow equations for n, K, and y. Equation 
( [4.51| ) for is identical to the equation for Kfj on a flat membrane p|,|l3| with K replace 



by K. Flow equations for K{1) and y{l) can be obtained in the usual way ||3^ from Eq. 
( [4.51| ). By breaking up the integral dr into J^^ dr + J^t dr, putting the first part into a 
renormalized K{1) and rescaling the second part by r ^ re~', we obtain 

iPKn)-' = mil))-' + —y\l) / drr'-^-P^ (5.1) 



p2 



with 



d{(3K)-' At 



3 



dl 



dy ( nPK v 



Observe that the bending rigidity does not appear explicitly in these equations; it only 
appears via the dependence of PK on k. This suggests that there should be an effective 
theory for a flat membrane with a bare hexatic rigidity equal to K = K — {3/327r){K/ k,)'^ 
rather than K. Indeed, if we integrate our height fluctuations to produce an effective flat 
space Hamiltonian via 



J Vhe-f'^--f'>^o, (5.4) 



we obtain 
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?ieff = IPK J d^xiVey. (5.5) 

Thus, the fluctuating membrane behaves hke a flat membrane with a hexatic rigidity K. 
Note that the renormahzed rigidity K appears in Eq. ( ^.5|) because fluctuations at all 
wavenumbers q have been integrated out. If we had calculated an effective Hamiltonian 
by removing h{q) for q within a shell A/b < q < A, there would have been no shift in K. 
The renormalization of k can be calculated using Eqs. ( |3.23|) and ( p.24| ). Diagrams 



contributing to the renormalized bending rigidity to lowest order in (3 ^ are shown in 
Fig. R|. They yield 
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_ _3_ /-^ ^ 3 -r^h^ {f^K'f ( dq ^ 

Air Jl-^ q Air 4(3k Jl-^ q Att 4/5k \Jl-^ q J 

_ 3 / (3K'\ r^dr 3 vrV (^^ / drV 

where L is the linear dimension of the membrane, which we take to be infinite, and A = 1/a. 
We can now follow exactly the same procedure to obtain the recursion relation for k we used 
to obtain those for K and y. To lowest order in (/3k) we can replace K' by K. We break 
the integral into two parts, and rescale to obtain 

dpK__3 (,_PK\ 

dl ~ AT^y A(3k)- ^^-^^ 

This equation was derived from the low-order diagrams shown in Fig. || in which the Coulomb 
energy K/ q^ is treated as a perturbation. One can also structure a perturbation theory in 
which the Coulomb propagator K' /q^ coupling S to S [Eq. (|3.24|) ] is replaced by K' /e^^ g^q^ 
where 

K' K 1 2K 

= 1 - -^Cssiq) - -^-^C^f^ + -jC^s (5.8) 

as shown in Fig. ^ Similar expressions can be derived for the effective Coulomb potentials 
coupling fi to n and n to S. In this scheme. Figs. §3 and c are replaced by Fig. |^. In 
addition, there are higher order diagrams such as those shown in Fig. These diagrams 
have two or more loops and are difficult to treat using our naive renormalization scheme. 
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In the companion paper, we will employ a controlled regularization scheme to treat these 
graphs using the sine-Gordon Hamiltonian. It shows that K should be replaced by K in Eq. 
(p.TI). Thus we have 

di ~ AT^y 4(3^) ■ ^^-^^ 

The renormalization flow equations now depend on K and not on K. 

Equations (|5.2|) , ( |5.3|) , and (|5.9|) define renormalization fiows for fiuctuating hexatic 
membranes. They have a fixed point at 

PK* = ^, y* = 0, p^* = p (5.10) 

We can obtain equations in the vicinity of this fixed point by defining 

I3K = l3K*{l-x), I3k = I3k*{1-z). (5.11) 



Then 



f^8.V, 1^2.. (5.12, 



dz 3 X — z 



dl 1f\-z 



(5.13) 



The fiow lines for these equations are shown in Figs. 10 and Fig. 11. Figure 10 a shows fiows 



in the a;|/-plane, which are identical to those in fiat space [ffl. Figure 10 b shows fiows in the 



yz-plane, which are similar to those in the xy-plane. Figure [TT|a shows fiows in the {^^K^~^- 
(/3k)~^ plane. As in previous treatments IQ], there is a fixed line in this plane at (5k = (3K/4 
{x = z), and the crinkled-to-crumpled transition occurs at (3K = j3K* independent of 13k, 
and the Kosterlitz-Thouless transition occurs at x = independent of z. 

The separatrix AO in Fig. p!0| a is the transition line dividing the crinkled hexatic phase 
from the melted crumpled phase. Systems whose initial values of x and y lie below and to 
the left of the separatrix fiow to the low-temperature crinkled hexatic phase; points above 
and to the right of the separatrix fiow away to the crumpled fiuid. The initial value of x is 
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Xo = 1 — K /Ti* , which increases with decreasing j3K = (3K — {3/327r){K/K)'^. Decreasing 
the bending rigidity k decreases K and eventually drives x to the right of the separatrix. 
Thus, there will be a transition from the crinkled to the crumpled phase as k is decreased. 
An alternative way to see this is to use K rather than K as the independent variable. To 
do this, we define K{1) to satisfy Eq. ([4.52|) with K and n replaced, respectively, by K{1) 
and k{1). Flows in the {{l3K{l))~^, {i3n{l))~^) plane are shown in Fig. |Tl|b. The vertical 
transition line in the {{l3K{l))~^ , plane is now curved and crosses the {(3K)~^ = 

axis at a finite value of {i3k)~^. Points to the left and below the separatrix flow to the 
crinkled fixed line. Points outside it flow to the crumpled phase. 

In deriving our recursions relations, we restricted ourselves to {i3K)~^ and 

(3K/{(3kY < 27r. Thus, our calculations strictly speaking only apply below the curves OA 
in Figs. ^ and |ll]b. The flows in Fig. |11] we calculated for p = 1. Figure ^ shows flows for 
the physically more interesting case of p = 6. Our approximation applies to a considerable 
region around the fixed point P. Higher order terms in {(3k)~^ and l3K/{f3nY would have to 
be included to get an accurate picture of what happens above the curve OA. Is seems likely 
to us, however, that these higher order terms will not lead to any qualitative modifications to 
Fig. |[ For any finite value of K"^ and k"^, the crinkled-to-crumpled transition is governed 
by the fixed point P. For infinite hexatic rigidity, = 0, the transition is controlled by 



different physics, probably analogous to that of tethered membranes |T0|. There could, of 
course, be some phase boundary at large f3K/{(3nY separating K = oo-like behavior from 
finite- ii" behavior, but we do not see any particular reason why this should be so. 
The persistence length in the fluid phase is 

^P = ae'\ (5.14) 

where I* is determined by 

(3k{1*) = (3k*{1- z{l*)) = 0, (5.15) 

i.e., by z{l*) = 1. Integration of Eqs. ( |5.12| ) and ( |5.13| ) yield z{l) and thus ^p. Equation 
( ^.13| ) for 2; is a nonlinear and cannot be determined analytically. We can, however, obtain a 
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very good estimate of when > ^kt, where ^kt ~ «exp(fe/|T— Tx^p/^) is the Kosterhtz- 
Thouless correlation length. The latter length is C,kt = ae^^ ■, where (3K{Iq) = or a;(Zo) = 1- 
The recursion relation, Eq. |5.9| for z{l) is only valid for K > {]. When 7^ — > 0, we approach a 
fluid phase with recursion relations for k determine by n alone. We, therefore, assume that 
x{l) = for / > Iq. Then 

dz ?) , , , , 

^ = v- ' > 

This equation can be integrated subject to the boundary condition that z{Iq) be the value 
of z{l) aX I = Iq determined by Eqs. |5]^ and |5.11| . Thus, 



z{l) = z{lo) + ^{l-lo). (5.17) 

/* = /o + 2p2(l-2(Zo))/3, and 

= ae'* = exTe^'^^''('")/^ (5.18) 

where we used (3k{Iq) = p^(l — z(/o)/(27r). This expression is to be compared with the result 
a pure fluid membrane. It is what one would naively have expected. The 
microscopic length a is replaced by the KT coherence length C,kt, and k at length scale a is 
replace by k at length scale ^kt- Equation ^.18| for is valid provided C,p > ^kt (or I* > Iq). 



Near the KT critical point, this inequality may not be satisfied. 



VI. REVIEW AND DISCUSSION 

We have investigated the Kosterlitz-Thouless transition on fluctuating hexatic mem- 
branes. We developed three equivalent Hamiltonians for describing these membranes: the 
hexatic elastic Hamiltonian expressed in terms of gradients of the hexatic angle variable, a 
Coulomb-gas Hamiltonian, and a sine-Gordon Hamiltonian. The Coulomb gas is character- 
ized by Gaussian-curvature and disclination charge densities, which interact via potentials 
partially determined by the Liouville action arising from a covariant cutoff. The sine-Gordon 
Hamiltonian has a linear coupling between the sine-Gordon field and Gaussian curvature. 
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We showed that height fluctuations, when integrated over all wave number, soften the hex- 
atic stiffness K. As a result, the disclination melting transition from the hexatic crinkled 
phase to the fluid crumpled phase is brought about both by increasing temperature and 
by decreasing the bending rigidity k. We derived renormalization group recursion relations 
for K, and the disclination fugacity y that explicitly verify this. In a companion pa- 
per, we provide an alternative derivation of these recursion relations using the sine-Gordon 
Hamiltonian. 

Though the picture we present of the Kosterlitz-Thouless transition from the hexatic 
crinkled phase to the crumpled fluid phase makes a great deal of sense, it does leave some 
incompletely answered or unanswered questions. First, we believe that the nature of shape 
fluctuation in the crinkled is not completely resolved. It is generally believed ||^ that the 
crinkled phase is characterized by power-law correlations in layer normals: (N(x) ■ N(0)) ~ 
|x|^'', where rj = 2T/{ttK). This result is obtained via exponentiation of the expansion 

(N(x) ■ N(0)) ^ 1 - |([V/i(x - V/i(0)]2) 

= l--^ln(|x|/a) (6.1) 
Ittk 

using K = K/4. The exponentiation of this series has not been explicitly justifled with, for 
example, a calculation of second order terms in {T/k) ln(x/a). It is interesting to note that 
the de Gennes-Taupin persistence length ^gt = ae^'^'^/'^ beyond which a fluid membrane 
is crumpled was calculated using Eq. (|0|) and setting (N(x) ■ N(0)) = 0. This observation 
would suggest that the possibility that the crinkled phase is in fact crumpled, but with a 
longer persistence length than the fluid phase, cannot be ruled out. 

Second, we have in our calculations imposed a constraint of charge neutrality. In flat 
membranes (i.e., fllms), this constraint is imposed by the prohibitive energy cost of having 
excess charge of either sign. We believe that there is a similar energy cost for breaking charge 
neutrality in free membranes, which are not constrained to be flat, though the situation here 
is more complicated. Nelson [^] has pointed out that the plus-minus symmetry present in 
flat membranes is broken in free membranes. Membranes with a single disclination undergo a 
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mechanical buckling transition from a flat configuration to a cone configuration for a positive 
disclination or a saddle configuration for a negative disclination. The energy of the cone is 
lower than that of the saddle, though both energies are proportional to ln(i?/a) where R is 
the linear dimension of the membrane. Does this asymmetry modify the picture presented in 
this paper? We believe not. A free membrane will choose configurations that will minimize 
its free energy. Any configuration, whether flat or buckled, with an excess of one sign 
of charge will have a contribution to its energy proportional to ln{R/a). Charge neutral 
configurations, on the other hand, have energies that are finite in the R ^ oo limit. Thus, 
there is a prohibitive energy cost in both flat and free membranes to the violation of charge 
neutrality. The shape of a membrane near the core of positive and negative disclination may 
nonetheless differ and lead to different fugacities and ?/_. We generalize our treatment 
to include this possibility in a companion paper. The results are that the ratio of the two 
fugacities is a marginal variable. In the ordered phase, both y+ and y_ scale to zero, and 
the KT transition to the disordered phase is not affected. In the disordered phase, both 
fugacities grow with a fixed ratio. A more complete treatment of Gaussian curvature will 
be needed to interpret this result. 

We are grateful to Mark Bowich, David Nelson, Phil Nelson, Burt Ovrut, and Leo Radz- 
ihovski for helpful conversations. This work was supported in part by the Penn Laboratory 
for Research in the Structure of Matter under NSF grant No. 91-20668. Further support was 
provided by the Institute for Theoretical Physics (under NSF grant No. 89-04035), where a 
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FIGURES 
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FIG. 1. Schematic phase diagram in the {PK)~^-{I3k)^^ plane adapted from Q. The vertical 
line BP is the Kosterlitz-Thouless disclination melting line of a flat membrane. The curved line 
PC is an estimate of the crumpling critical line based upon a comparison of the energy and 
entropy of a single positive disclination in a buckled membrane. The line OP is the mechanical 
buckling instability line of a membrane with a single positive disclination. Above this line the 
zero-temperature membrane is buckled. Decreasing K or increasing T at fixed k in the vicinity of 
P leads to a disclination melting transition to the crumpled phase. On decreasing k, at fixed K, 
the crumpled phase can only be reached by crossing the crumpling line. 
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FIG. 2. Phase diagram for hexatic membranes {p = 6) in the {(3K) ^-{(3k,) ^ plane obtained 
from Eqs. Q, Q, and Q. The dark curved hne [ (PK)-^ + {3/32tt){Pk)-^ = tx 112 = 0.0436] 



is the critical line separating the crinkled from the crumpled phase. The line OP is the crinkled line 
with 4(/3k)"^ = {PKy^ + (3/327r)(/3K)~2]. P is the crinkled-to-crumpled fixed point. Unlike Fig. 
1^, the critical line has curvature at P. The dashed curve OA is the curve 2-k{I3K)^^ = {f3n)^'^. The 
calculations in this paper are approximately valid in the region below this line where (/3k) < 1 and 
fiK < 2tt(I3k,)'^. The crinkled-to-crumpled transition occurs via disclination melting in the vicinity 
of P when either k or is decreased or when temperature is increased. Though our calculations 
do not apply above the curve OA, it plausible that the crumpled-to-crinkled transition occurs via 
disclination melting for all k and K except at K = oo, i.e., that the only effect of higher order 
terms in {(3k)~^ and f3K/{f3K)'^ is to change the shape of the critical line. Alternatively, there 
could be some other phase boundary above the curve OA separating melting from some kind of 
crumpling transition. 

H vAAA^ 
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(g) 

FIG. 3. Elements for constructing a diagrammtic perturbation expansion for a coulomb gas on 
a fluctuating membrane: (a) The height propagator Ghh = '^/{Pi^q'^), (b) the height propagator 
with two gradients represented by the vertical lines, (c) the bare vortex charge-density propagator 



C5- propagator [Eq. ( 4.44 )], (d) the coulomb vertex K' /q^ in TLss-, (e) the coulomb vertex —K/q^ 
in TLnSi (f) the coulomb vertex K/q^ in TLnn-, and (g) a nonlinear vertex from the curvature energy 
[Eq. dO)]. 
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(c) 



+ 





(e) (f) 

FIG. 4. Diagrammatic contributions to the polarization bubble V represented by the 
cross-hatched diagram to the left of the equal sign. The components of these diagrams are de- 
fined in Fig. 1. The leading order contribution to V is given by diagram (a). Diagram (e) is the 
first correction to V arising from the ^/g factors in the definition of V. The double crossline on the 
height propagators represent the gradients, with appropriate symmetry, in the Gaussian curvature 



[Eq. gl2D ]. 




FIG. 5. Diagrammatic representation of the full C?^ propagator when TinS = 0. 




FIG. 6. A representative additional diagram contibuting to V when UnS is turned on. The 
expansion for C^^ is still given by Eq. ( 4.44 ) and Fig. ^. 
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(b) 



(c) 

FIG. 7. Diagrammatic expansion of (a) Cnh, (b) (^^5, and (c) C-^. The bubble is C^^^ repre- 
sented in Fig. 0. 
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FIG. 8. Diagrams for the bending rigidity k 





(a) (b) 

FIG. 9. Diagrams contributing to k in terms of the screened Coulomb propagator (represented 
by the double dashed line) K'/{eu ssQ^) defined in Eq. (| 
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(b) 

FIG. 10. (a) Renormalization flows in the xy-plane. These are identical to the order of our 
calculations to those of the flat space xy-model. The separatrix AO is the critical line, (b) Flows 
in the zy-plane. They are similar to those in the xy-plane. 
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FIG. 11. (a) Renormalization flows in the xz plane. The transition from the crinkled to the 
crumpled phase takes place at x = independent of z (i.e., independent of j3n), in agreement with 
previous calculations Q. (b) Flows in the {l3K)~^-{f3n)~'^ plane. Here, there is a transition from 
the crinkled to the crumpled phase as k is decreased. 
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